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^ ■ ABSTRACT 

o 

This paper investigates the dynamics of excitonic transport in photocomplex LHCII, the primary component of 
the photosynthetic apparatus in green plants. The dynamics exhibits a strong interplay between coherent processes 
mediated by the excitonic Hamiltonian, and incoherent processes due to interactions with the environment. The 
. spreading of the exciton over a single monomer is well described by a proper measure of derealization that allows 
one to identify two relevant time scales. An exciton initially localized in one chromophore first spreads coherently to 
neighboring chromophores. During this initial coherent spreading, quantum effects such as entanglement play a role. 
As the effects of a decohering environment come into play, coherence and decoherence interact to give rise to efficient 
and robust excitonic transport, reaching a maximum efficiency at the levels of decoherence found in physiological 

Cln' conditions. We analyze the efficiency for different possible topologies (monomer, dimer, trimer, tetramer) and show 

+Jj . how the trimer has a particular role both in the antenna and the wire configuration. 
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INTRODUCTION 



Recent experiments probing excitonic transport in green sulphur bacteria suggest that quantum coherence plays 
an important role in photosynthesis 0-0] - Detailed models of the interplay between coherent exciton dynamics 
and decoherence and relaxation induced by the exciton's environment show that the resulting transport is robust 
and efficient [EH?}, an effect known as environmentally assisted quantum transport (ENAQT). This paper extends 
these analyses to the light- harvesting complex of green plants [Hf , specifically, transport through the light-harvesting 
complex LHCII . Excitonic transport through sets of coupled LHCII complexes |8l4ll] differs in significant ways from 
• the transport through the Fenna-Matthews-Olson complex (FMO) of green bacteria . Notably, the LHCII can act 
both as antennae and 'wires' capable of transferring the excitons captured by other complexes through the structure. 
In order to move through a sequence of LHCII complexes, the exciton must move both up and down in energy, a 
. . , process mediated by interactions with the environment. LHCII complexes can be found in different forms allowing for 
various regulation activities. In particular, some LHCII complexes can migrate, under specific light conditions, from 
the Photosystem II (PSII) to Photosystem I (PSI) in order to optimize the photosynthetic process (state transition), 
or disassamble into monomeric subunits in order to favor the regulation of light harvesting in excess light [14j . 

This paper applies a purely-decohering Haken-Strobl model to analyze interactions between coherent and decoherent 
dynamics in excitonic transport in LHCII 12, T^. The advantage of using Haken-Strobl model is that it is the simplest 



model that allows this transport to be investigated in the strongly-coupled, non-perturbative regime. The disadvantage 
of a purely decohering model is that it does not include relaxation, and so will over-estimate transport rates as the 
exciton moves up in energy, and under-estimate them as it moves down. Nonetheless, our analysis shows that pure 
decoherence is a surprisingly effective transport mechanism even when the exciton is moving down through an energy 
funnel. 

To analyze the interplay of coherence and decoherence, we use familiar tools from quantum information theory. 
Quantum mutual information between sets of sites is used to track the pathways by which correlations spread through 
the LHCII monomeric complex due to exciton motion. Negativity and concurrence are used to demonstrate and 
quantify entanglement. The comparison among the different measures allows us to identify two timescales. Over 
the first half picosecond, an initially localized exciton spreads coherently to neighboring chromophore. This coherent 
spreading is accompanied by rapid oscillations in quantum mutual information and negativity, indicating the presence 
of significant coherence and entanglement. Next, over the course of several picoseconds, the coherent oscillations 
disappear, the negativity decreases, and quantum mutual information between sites grows as decoherence kicks in and 



FIG. 1: The chains of Chls in the trimeric LHCII. Blue: CHls-b; green: CHls-a. 



the exciton diffuses partially incoherently to more distant sites. The interplay between coherence and decoherence gives 
rise to highly efficient transport through the complex. We analyze the efficiency of transport for various geometries of 
LCHII (monomer, dimer, trimer, and tetramer) and for the two possible configurations: antenna and wire. We show 
that efficiency increases as the number of subunits increases, saturating at the trimer. That is, the dimer exhibits 
more efficient transport than the monomer; the trimer is more efficient than the dimer; and the trimer and tetramer 
give the same efficiency. 



BACKGROUND 



Photoabsorption, the initial step of photosynthesis, takes place in photosynthetic complexes formed by groups of 
pigments (chlorophylls, Chi) and proteins placed within the thylakoid membrane The light-harvesting pigments 
are arranged in protein matrices in such a way that the photo-excitation is funneled to the reaction center, where the 
energy carried by the excitation is used for fueling chemical reactions. 

One of the main elements in the photosynthetic complexes of higher plants is the light-harvesting complex LHCII. 
In the following we will refer to the LHCII that can be found in the Arabidopsis thaliana, which is a flowering plant 
that has been extensively studied as a model in plant biology and genomics. In this plants, LHCII can be found as 
peripheral antenna in the supercomplex PSII 13]. Its role in the supercomplex is to collect photons from the incoming 
light or excitons from the neighbouring complexes (other LHCIIs or other minor antenna complexes) and transfer 
them to other LHCII complexes or to the reaction center. 

The LHCII crystal structure has been determined up to 2.5 A in resolution lg, [l7|: it is composed of three similar 
monomeric subunits (Lhcbl-3) each containing 14 Chls molecules embedded in a protein matrix (see FigfTJ). Two 
different types of Chls are present: 8 Chl-a type, and 6 Chl-b type. The main difference between the two is the 
presence in Chl-b of a carbonyl group that allows for higher excitation energies. The Chls are disposed in two layers 
which are called stromal and lumenal, the first being oriented towards the outer part of the thylakoid membrane, the 
second towards its inner part (see fig. [2]). 

Within each layer, the Chls can be grouped on the basis of their relative distance and orientations which determine 
the strength of the electronic interaction. The main groups of Chls are depicted in fig. |3j different groups have 
different linear absorption spectra and two main bands can be identified: the Chl-a band centered at 14925 cm" 1 , and 
the Chl-b band centered at 15385 cm" 1 . 

The exciton energies and the relative pigment participations in LHCII have been determined in The energy 
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FIG. 2: LHCII monomer in the membrane, and its top and bottom view. Blue: CHls-b; green: CHls-a 



relaxation pathways of the system were described by van Grondelle (see [8J) and refined in [10| by means of the 
experimental data obtained by applying 2D spectroscopic techniques to samples of LHCII at 77K; in the same paper, 
the Hamiltonian for the monomeric unit we will use in the following was also derived and optimized in order to 
account for the experimental results. The Hamiltonian is written in the site basis: the diagonal terms are the site 
energies, determined by fitting the linear absorptions spectrum and the off-diagonal terms account for the Coulomb 
interactions between pairs of Chls (see [lOj and references therein). The experimental evidences allows to identify two 
main downhill relaxation pathways [H, llCf . One takes entirely place at the stromal level (stromal-stromal) while the 
second starts in the lumenal level (lumenal-stromal) (see fig. U]). 

One of the fundamental mechanisms that is at the basis of the energy transfer through the LHCII is given by 
the interplay between strong electronic couplings between nearby Chls. The coupling allows both for the energy 
splitting between excitonic levels and for the exciton derealization. If we focus for example on the group of Chls 
a610, a611, a612 (fig. [3]) we see that the contiguity and the relative orientation of these molecules result in a strong 
interaction which in turn it allows for the presence of three exciton levels (levels 5, 2, 1 in fig. 2]) which are separated in 
energy and at the same time are spatially overlapping. In presence of exciton-exciton coupling, for example mediated 
by the environment, the spatial overlap allows for fast (< 100/s) relaxation processes within the group-excitonic band 
(levels 5,2,1 in fig. gl). 

The same mechanism can be seen in another group of Chls 6601,6608,6609 which allows for a fast (< 100/s) 
relaxation process within the stromal b-band (levels 13, 11, 10 in fig. @|. 

In analogy with what happens for the previous groups of Chls, the intra-group energy pathways can be described 
by the same mechanism. In particular, in the lumenal-stromal pathways, one can identify a bottleneck of the pathway 
that is given by the Chls 6605 and a604; here the interaction between the two molecules gives rise to a pair of well 
localized excitonic states (levels 8 and 9 in fig. 0|) which have poor spatial overlap with the lower energy states 
concentrated on the Chl-a groups in the lumenal and the stromal side. The localization of the wave function gives 
rise to the experimentally determined slowdown in the relaxation process relative to the lumenal-stromal pathway. 
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FIG. 3: Stromal (A) and lumenal (B) layers of the LHCII. The groups of strongly coupled Chls are enclosed in thick circles. 
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FIG. 4: Energy relaxation pathways as experimentally determined in [1C 



The dynamics of the bottleneck will be seen in action in the simulations of exciton transport through the monomer 
in the next section. 

The lowest energy states of the monomer are two of the excitonic states localized in the et610, a611, a612 group of 
Chls. These sites are therefore indicated as the output sites of the monomer unit (donors), which can be coupled 
to other LHCII monomeric units, or other complexes in the PSII. The coupling can be between exitonic states 
corresponding to similar energies but located on different neighbouring complexes. In (lQj these two output states have 
been studied in terms of their directionality, which is determined by the interplay between the site basis contribution 
to the specific excitonic state (derealization over the group of Chls a610, a611, a612) and the relative orientation of 
the donor transition dipoles. The result of the calculations showed that the two excitons can act as donor states in 
two different directions, and this mechanism has been suggested as a way that LHCII complexes have to optimize the 
exciton transfer to other complexes even in presence of their misalignment. In our analysis we want to compare two 
different ways in which the monomer can be coupled with the outer complexes: i) a sites based coupling, where each 
output site (Chls alO, a611, a612) is independently coupled with an external complex that we model as a sink; i%) an 
exciton based coupling, where the two lowest excitonic states are independently coupled with an outer sink. 

For the inter-monomeric couplings, the experimental evidences in (l7j (fioj). show that the main coupling should 
be localized in the stromal side between Chls of b type (&601, 6608, &609); the strength of the coupling should be of 
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the order of 42cm -1 , (35cto _1 ). In [Hj this inter- monomeric coupling has been neglected in the determination of the 
single monomer Hamiltonian. While this has the effect of shifting the excitonic b-band, it should not significantly 
affect the other transitions which are localized within each monomer. In order to have a description of the whole 
trimeric LHCII, we reintroduce the coupling between Chl-b pertaining to different adjacent monomers; in particular 
we choose two different configurations: 6601 <->■ 6609 and 6601 f-> 6609, 6608. 

We finally note that another Hamiltonian for the LHCII complex has been derived in [ll|. However the results that 
can be obtained in the following analysis by using the data reported in [ll| instead of the ones reported in [l(J do not 
change in a significant way. 



MODEL AND TOOLS 



In the following we will focus on transport properties of the monomeric, dimeric, trimeric and tetrameric LHCII 
systems. We study the dynamics of the open quantum system in the presence of dephasing processes (Haken- Strobl 
formalism) in the presence of recombination and trapping mechanisms. In this model we resort to the description used 
in @, @] , where recombination and trapping processes are modeled by adding non-Hermitian terms to the Hamiltonian. 
The effects of static disorder will be taken into account below. The equation of motion for the density matrix of the 
monomer subunit can be written as 

~tT [-^ ^ monomer j ^ J dephP(^} {Elrecomb ~l~ -^trapping ) P\ — 

' = L[p(t)}, (1) 

The free Hamiltonian of the monomer is a tight binding Hamiltonian and it is expressed in terms of the site energies 
and couplings given in (Io| 

onorner 

= Y^ e rn\rn)(m\+^2v mn (\m)(n\ + h.c) (2) 

m m<n 

The term 

L dephP (t) = j^AmPi^Al - ~ {A m Al,p(t)} , (3) 

m 

accounts for the presence of pure dephasing, A m = \m)(m\. The term 

Hrecomb = -tT ^ \m) (m\ , 

ni 

accounts for the recombination processes. For the trapping, we suppose that once the exciton has reached the output 
sites a6I0, a611, a6I2, it leaves the LHCII complex with a rate kt ra p- The trapping process can be expressed either 
with respect to the site basis i.e., the sites a610, a611, a612 are supposed to be singularly linked with other surrounding 
complexes 

Htrapping = iktrap ^ ' (4) 

m=610,611,612 

or with respect to the two lowest exciton eigenstates that are localized in the a6I0, a611, a6I2 sites: 

Htrapping = —iktrap(\El){Ei \ + \ El )(-E?2 | ) • (5) 

The main difference between the two pictures should be that with the site-based trapping mechanism the sink acts 
on the output sites, including the case where only the highest exciton localized on the output states (number 5 in 
figure 2]) is populated. By contrast, with the excitonic-based mechanism only the two lowest energy excitonic states 
are involved in the trapping dynamics. 

The functionals we use to evaluate the efficiency of the transport are the efficiency 77, defined as 



i] = 2kt rap y ^ I 

™ Jo 



dt{m\p(t)\m) (6) 
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FIG. 5: Energy levels for the Hamiltonian of the monomer and the trimer. 



and the average transfer time r, defined as 



kfrap ^ ^ I 

™ Jo 



dtt (m\p(t)\m) 



(7) 



The system dynamics and the efficiency 77 can be evaluated numerically by vectorizing the density matrix and 
constructing the proper linear super operator associated to L 

P(t) -> \p(t))) 



dp® 
1 dt 



)) = CW))) 



p(o))). 



C has been constructed using the identity \ABC)) = A <E> C l \B)), where A,B and C are matrices of size n and |-)) is 
a vector of size n 2 . We can now compute 77 in terms of C as 
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^fc^HzrXo))). 



Before passing to a naly ze the results of our simulations a comment on the monomer Hamiltonian is in order. The 
Hamiltonian given in [10( has been obtained after many optimization processes with the goal of faithfully reproducing 
the 2-D spectroscopy experimental results. This implies that the actual eigenstates that can be derived by diagonalizing 
the Hamiltonian are delocalized over groups of Chls which are sometimes different by the ones showed in fig. |4j where 
the relative pigment participations are derived via other optimization processes. The differences manly involve the 
highest energy eigenstates of the b lumenal and stromal band. In our analysis, we stick to the experimentally optimized 
Hamiltonian and use its actual eigenstates; the spectrum of the monomer and the trimer is given in figure [5] The 
maximal, minimal and average difference in energy for the monomer are A m i n w 4cm _1 , A max rs 172cm~ 1 , Aa v — 
59cm" 1 (ksT w 200cm~ 1 at room temperature). The differences in energies for the dimer and the trimer are very 
similar. 

In order to study the dynamics of the monomer we make use of well established quantum information measures of 
correlations. In particular, we measure the total amount of correlations between two subsystems of chromophores A 
and B with the quantum mutual information [201 ] 
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Tab = (S)a + (S) B - (S)ab (8) 

where (S)x = X); 1°§2 Aj is the von Neumann entropy of the reduced density matrix px of subsystem X evaluated 
in terms of its eigenvalues {A^}. As for the quantum correlations between subsystems composed by arbitrary number 
of chromophores we use the negativity Nab pl| that in the single exciton manifold can be written as 22] 



fc N 

a oo + 4 Y I ( n \PAB\m) | 2 - a 00 (9) 

n=l ra— fc+1 

where aoo is the element corresponding to the zero exciton subspace and A — 1, • • • , fc, B = fc + 1,-- - ,N are two 
generic subsystems of chromophores. The quantum correlations between pairs of sites m,n are also measured by 
means of the concurrence [lH which in the single site exciton has the simplified form 

Cm,n — 2|p mTn |. (10) 

In order to study the relationship between the dynamics of the above described correlations measures and the 
derealization of the exciton over the monomeric structure we define a measure of derealization T>(t) that involves 
the single site populations. In keeping with this paper's methods of using information-based measures to characterize 
the excitonic transport, we use Shannon entropy as a measure of derealization. If (m) is the population of site i at 
time t, then by using the normalized populations A^ = (n*) / ^2 ■ (n*) we can define 

V(t)=-Y,^og\i (11) 

i 

as the Shannon entropy of the normalized populations. The higher T>(t), the flatter the probability distribution {A^}, 
the higher the derealization of the exciton over the complex. 
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MONOMER DYNAMICS 



In the following we describe the dynamics of the monomer. The values of the recombination and trapping coefficients 
used for the simulations are the ones used for the FMO complex in The recombination coefficient T = 10~ 3 ps _1 

takes into account the estimated lifetime of the exciton, Ins. The trapping coefficient is k tra p = lps -1 and is assumed 
to be equal for each output exciton state. Our results do not depend sensitively on the exact value of the exciton 
lifetime and the trapping rate: what is important for the analysis is that the exciton has a relatively long lifetime 
compared with coupling rates, and that the trapping rate is comparable to those rates. 

We first focus on the time simulation of the evolution of the monomer in order to identify the possible energy transfer 
pathways. We therefore fix the value of the dephasing rate 7^ « 3ps _1 that corresponds to 77 K (temperature at 
which the experiments were done in [10|). The dephasing rate can be written in terms of the bath correlation function 
as 0: 

l(j>{T) = 2n^fd u J(u>)\ a = = 2tt^^- (12) 

where we have supposed to have an Ohmic correlation function J(w) = 7r fcB ^- Er j*- expw/u; c (super- and sub-Ohmic 
correlation functions give a similar dependence on T, E r , and lu c ). The recombination energy E r = 35cm -1 and the 
cut-off frequency uj c = 150cm -1 are chosen to be the ones used for the FMO simulations. Again, the qualitative 
behavior of the excitonic transport does not depend sensitively on the precise values of E r and uj c . 



Energy and correlation pathways 

We want to describe the time evolution of the state of the monomer coupled with the environment. We first try to 
identify and characterize the existence of two possible pathways, stromal-stromal and lumenal-stromal, by which the 
exciton, starting from a high energy state belonging to the b band, can reach the output sites. We therefore focus our 
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FIG. 6: Schematic picture of the bipartitions used for the analysis of correlations dynamics 

attention on the behaviour of the system without any trapping and for a value of the dephasing rate that corresponds 
to the temperature of 77 K used in [Io| . As noted above, since the Haken-Strobl model is purely dephasing and does 
not include an explicit relaxation term, we expect our analysis to underestimate the rate of transfer from high energy 
states to low energy states. Nonetheless, as will now be seen, pure decoherence gives efficient excitonic transport 
down the energy ladder. 

The distinct pathways can be studied by plotting the following populations: i) the populations of the excitonic 
states which are mostly localized on Chl-b molecules that belong to the stromal (lumenal) side Pbstrom{PbLum) H) the 
populations of the excitons states which are mostly localized on Chl-a molecules that belong to the stromal (lumenal) 

side PaStrom{,PaLum) • 

The different excitonic behaviors along the two distinct pathways can also be highlighted by using measures of 
quantum correlations. In particular, we study the quantum mutual information between the relevant subsystems that 
are naturally suggested by the energy landscape in fig. 01 in this way we can identify the correlation pathways and 
their dynamics. On the lumenal side we select the subsystems bL = {606, 607}, abL = {605, 604}, aL = {613,614}, 
while on the stromal side we select the subsystems bS — {601, 608, 609}, aintS = {602, 603}, aoutS = {610, 611, 612} 
and aS = aintS U aoutS; the bipartitions are schematically depicted in fig. [6] The growth of quantum correlations 
between subsystems is a signature of the spreading of the initially localized exciton between subsystems. The form 
that these correlations take over time reveals the mechanism of this spreading - an almost purely coherent initial 
propagation followed by semi-coherent diffusion. 

For the stromal-stromal pathway we choose as initial state of our simulations the highest energy eigenstate I-E14) 
of the Hamiltonian, which is mostly localized on the 6601 Chi. In Fig|7]we see how the exciton mostly flows from 
the stromal b-band to the stromal a-band on a very short time scale (~ 5ps). On a slower time scale the population 
partially delocalizes over the lumenal band. The flow of population between the two b-bands was highlighted in the 
energy pathway given in Q (but not highlighted in [Io|). where the possibility of a flow from the b-lumenal to the 
b-stromal band is estimated to be of the order of 2ps. Here we observe the inverse passage b- stromal to b- lumenal 
and this is due also to the partial derealization of the \E\a) on the b-lumenal sites. The global population decreases 
because of excitonic decay with a time scale of the order of Ins. 

The dynamics can be further analyzed by considering the correlation pathways. Figure [5] refers to the stromal- 
stromal pathway starting with I-E14). The plots shows that the dynamics mostly takes place on the stromal side. In 
particular (left plot) the bS subsystem initially gets correlated with a-band stromal sites as a whole (MlaSbS); in the 
first few picoseconds most of the correlations are established between the subsystems bS-aintS, and the subsystems 
aintS-aoutS. From the first picosecond on the bS sites get directly correlated with the output sites aoutS. The 
right plot in fig. |8] shows that while there are correlations between the stromal and the lumenal b-bands (MlbLbS), 
the correlations among the subsystems on the lumenal side and the lumenal-stromal correlations at the level of the 
a-bands are negligible (~ one order of magnitude smaller). This picture is consistent with a dynamics mostly localized 
on the stromal side of the monomer. 
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FIG. 7: Stromal- Stromal pathway: populations for initial state \Eia) 
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FIG. 8: Stromal- Stromal pathway: mutual information for initial state |i?i4); left: stromal side; right: lumenal side 



We now pass to analyze the behaviour of the system when the exciton starts on a a high energy eigenstate I-E13) 
which is mostly localized on the lumenal side (in particular on the site 6606). The behaviour of the populations shown 
in fig. [9] accounts for the presence of a bottleneck in the energy pathways @, EJ. The latter is due on one hand to 
the high spatial localization of the exciton involving the site a604 (which is far from the other a- lumenal sites), and 
on the other hand to the large energy separation with respect to the excitons localized on the neighbouring lumenal 
sites 6605,6606,6607 (AE « 172cm -1 eigenvalue 8 with respect to eigenvalue 9 in fig. [5]). Indeed, while in the first 
few picoseconds, the population of the b-lumenal sites decreases in favor of the population of a-lumenal band, part of 
the b-lumenal population start to flow toward the b-stromal band. The overall effect is that the a-stromal sites, and 
therefore the output sites, become populated with a smaller rate than in the stromal-stromal case. 

The trapping effect on the lumenal side can also interpreted in terms of correlation pathways. As it is shown in the 
right plot of fig. [10] the b-lumenal band is initially highly correlated with the b-stromal band (black line, right plot). 

On one hand this is due to the fact that the eigenstate \Ei3) is partially delocalized on the sites 6608,6609. On 
the other hand, as already mentioned, the b-lumcnal to b-stromal flow, was already pointed out in terms of energy 
pathways in [8[ , where the transfer time between b-lumenal sites and b-stromal sites were estimated of the order of 
2 — 4ps, and, despite the small level of interaction (< 11cm -1 , see Ham) it is consistent with the contiguity of the 
b-lumenal and b-stromal sites in the monomer. 

While the correlations between b-lumenal and b-stromal bands rapidly decay in the first few picoseconds, there is 
not a corresponding growth of correlations on the lumenal side: the b-lumenal and a-lumenal sites remain very poorly 
correlated among each other and with the rest of the a-stromal sites. At the same time there is a clear enhancement 
of the correlations on the stromal side, which become rapidly greater than those on the lumenal side, suggesting the 



10 



0.8 




5 10 15 20 25 30 35 40 45 50 

ps 



FIG. 9: Lumenal- Stromal pathway: populations for initial state I-E13) 




FIG. 10: Lumenal-Stromal pathway: mutual information for initial state IS13}; left: stromal side; right: lumenal side 



activation of the stromal-stromal pathway. 

The above analysis is robust with respect to the choice of the initial state; for example in fig. [TT]the same simulations 
have been carried out for an initial state |6607) fully localized on the b-lumenal site 6607. Here the activation of the 
stromal-stromal pathway is evident; in the first lps the b-lumenal subsystem gets correlated both with the sites 
6605, a604 (M I ab Lb L) and with the b-stromal band {MlbLbS right plot); then the correlations are mostly built on 
the stromal side (left plot) while on the same time scales the correlations with the output sites on the lumenal side 
are build with a slower rate. 

The net effect of the presence of the above mentioned bottleneck on the lumenal side is a substantial slowdown 
of the lumenal-stromal dynamics. There are a number of possible functional advantages for this slowdown. One 
possibility is to assist in photoprotection, the elimination of triplet exciton states that can create harmful singlet 
oxygen. This elimination takes place primarily by the transfer of triplet excitons to triplet carotenoid states, and has 
a relatively slow timescale (a fraction of a microsecond) compared with singlet exciton transfer [l9| . Triplet transfer 
is a Dexter process, mediated by wave function overlap, and requires the carotenoids to be physically close to the 
chromophore carrying the triplet. X-ray crystallography studies of LHCII suggest that the a604 chromophore in the 
lumenal bottleneck is close enough to a lutein carotenoid to allow triplet transfer, a process confirmed by spectroscopy 



A second possible function for the lumenal bottleneck is that the bottleneck 6605 chromophore could mediate 
excitonic transport from one trimer to another [t7| . This chromophore 'sticks out' from the others in the LHCII 
crystallographic structure, giving both weaker couplings to the other chromophores within the LHCII monomer, 
and potentially stronger couplings to chromophores in neighboring triples. Crystallographic investigations of LHCII 
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FIG. 11: Lumenal-Stromal pathway: mutual information for initial state |6607); left: stromal side; right: lumenal side 



suggest that the 6605 is positioned to transfer excitons to the 6606 and a614 chromophores of neighboring trimcrs 



in LHCII aggregates jl7| . This transfer pathway could also participate in photoprotection via non-photochemical 
fluorescence quenching. 

A third possible reason for the slower lumenal pathway is that it might allow two excitons to propagate through 
the LHCII complex simultaneously without quenching. The weak coupling between the bottleneck chromophores of 
the lumenal pathways and the chromophores of the stromal pathway, together with their relative spatial separation, 
could allow an exciton localized in the lumenal pathway to wait for a stromal exciton in the same complex to pass 
through the stromal trapping states, before passing through itself. 

To summarize, the dynamics of the slow lumenal side and the fast stromal side have many have a rich set of 
potential biological functions. These dynamics are in turn based on a rich quantum structure, which the next section 
elucidates. 



Delocalization and quantum correlations 

In the following we examine the monomer dynamics from the point of view of the time scales that characterize the 
delocalization of the excitons over the whole monomeric structure and the quantum correlations between the various 
subsystems. In order to estimate the delocalization time we plot the populations of the various bands in fig. [12] 
(left plot), and the previously denned measure of delocalization T>(i) (fTTj) . Since we want to study the dynamics of 
the spreading of the correlations we first focus on the initial state |6601) localized on the Chi 6601 only. We choose 
this state because it is very close to the highest energy eigenstate -E14}, which is mostly localized on the same site 
but has non-negligible quantum correlations with the rest of the system. We want to start with a state localized on 
a single chromophore, in order to study how the correlations spread through the structure. In fig. [12] we plot the 
delocalization function V(t) for 7^ = 3ps 77K), T = 0.001 and no trapping. The delocalization has a very fast 
growth and it can be well represented by a function V(t) = yO + Al exp — t/t\ + A2 exp —t/t? where two time scales 
t\ » 250 fs and ti ~ 2.68ps appear. The relevance of these time scales can be understood by studying the quantum 
correlations between the relevant subsystems. 

In the left panel of fig. [13] we plot the mutual information for the relevant subsystems over a time of lOps. As it 
can be easily seen, the second time scale £2 ~ 2.68ps can be correlated with the growth of the mutual information 
between the various subsystems which reach its maximum at a time t very close to ti. In particular the correlations 
between the b-band and the a-band on the stromal side has a maximum at t 2.63ps. In the right panel of fig. I13l we 
show the dynamics of correlations within the first picosecond. Here the dynamics displays an initial fast growth of all 
correlations (« 0.50ps) which are characterized by an oscillating behaviour; the time scale of these oscillations is of 
the order of O.lps. These oscillations are signatures of the high degree of quantum coherence in the initial spreading. 
The growth becomes regular within the first 400/s ~ 1.5ii; within the same period of time the correlations spread 
toward the lumenal side (MlbLbS) 
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FIG. 12: Initial state |6601) (« \Eu)), state localized on the Chi 6601 on the stromal side. Populations (left) and derealization 
(right) T>(t) with fitting curve. 
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FIG. 13: Initial state [6601) (~ \Eu)), state localized on the Chi 6601 on the stromal side. Mutual information between 
subsystems on the stromal side and between b-lumenal and b-stromal band (MlbLbS); the right plot displays the dynamics 
during the first lps. 



The overall effect can therefore be described in terms of the two relevant timcscales: the derealization takes place 
over the structure with an initial fast rate, and subsequently it grows toward its asymptotic behaviour with a slower 
rate. The time scale of the initial rapid and oscillating growth is consistent with the dynamics of the quantum 
correlations present in the system. In fig. [14] (left plot) we show the behaviour of the negativities among the relevant 
subsystems over the first lOps (NegbLbS is the negativity between the lumenal and the stromal b-bands). In the first 
picosecond (right plot), after a first rapid growth (« 50/s ) they show the same oscillatory behaviour of the mutual 
information and they then decay in a smooth way after the first « 400ps (left plot). This behaviour is also shown by 
the concurrences between sites: in fig. [T5]we show the relevant (non-negligible) concurrences between the site 6601 
and other sites. Non-zero negativity and concurrence demonstrate the presence of entanglement during the initial 
coherent spreading, similar to the presence of entanglement in FMO 22|, |24j. The quantum correlations are therefore 
established in the first few hundreds of femtoseconds and they give their contribution for the first rapid growth of the 
derealization of the exciton over the whole structure (stromal and lumenal side). 
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FIG. 14: Initial state |6601) (~ \Ei4)), state localized on the Chi 6601 on the stromal side. (Un-normalized) Negativities 
between subsystems on the stromal side and between b-lumenal and b-stromal band (NegbLbS); the right plot displays the 
dynamics during the first lps. 
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FIG. 15: Initial state |6601) (~ |_Ei4)), state localized on the Chi 6601 on the stromal side. Relevant, i.e. non-negligible, 
concurrences between site 6601 and other sites on the stromal (6603,611) and lumenal (6606) sides; the plot displays the 
dynamics during the first lps. 



Effect of the trapping on the monomer dynamics 

In order to see how the trapping modifies the monomer dynamics we focus on the stromal-stromal dynamics. The 
initial exciton state is I-B14) and we choose to fix the value of the dephasing rate to 70 = 12ps~ 1 (that approximately 
corresponds to the ambient temperature if one uses [l~2|) . In fig. [TH the left plot displays the populations of the 
various subsystems. The comparison of the populations dynamics with the site based trapping mechanism and the 
exciton based mechanism shows that the first one is obviously more efficient in reducing the total population and the 
population of the various band (in particular the a-lumenal one). As noted above, this difference arises in our model 
because site-based trapping operates on three sites, while the exciton based mechanism acts only when the two lowest 
eigenstates get populated. 

In the right plot of fig. Ql)]we show the derealization functional for no trapping and for site/exciton based trapping. 
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FIG. 16: Initial state |2?i4). Dephasing 7^ = 12ps _1 ; trapping kt rap = lps -1 ; recombination T — O.OOlps -1 . Left: populations 
of the different bands with exciton/site base trapping; right: delocalization functional T>(t) with exciton/site base trapping and 
no trapping 

The presence of the trapping mechanism starts to be relevant already after the first picosecond, i.e. before the full 
delocalization of the exciton has occurred. Indeed for 70 = 12ps _1 the typical time scales of the delocalization without 
trapping are t\ « 400/s and ti ~ 1.5ps. This means that, in presence of trapping, the delocalization due to the initial 
dynamics of the quantum correlations has a fundamental role in the energy transfer process. 

The monomer as wire 

An interesting problem is to determine the typical time scales that govern the delocalization of an exciton initially 
localized on the output sites a610, a611, a612 of the monomer. These quantities become relevant when one describes 
the behaviour of the dimer and the trimer when they act as quantum wires. Indeed, the LHCII complex can in 
principle be activated by other neighbouring LHCII complexes and this should happen when the output sites of a 
pair of LHCII arc sufficiently close. In this process, the a610, a611, and a612 sites on an LHCII monomer accept 
an exciton from the same sites on a neighboring trimer. The exciton then spreads first through the monomer, and 
then throughout the three LHCII units of the trimer. When it reaches another set of output sites, the exciton can be 
transferred to another trimer, and the process repeats. 

Accordingly, we now analyze the interaction between donors a610, a611, a612 on one monomer within the trimer, 
and "acceptor" a610, a611, a612 sites on a second monomer within the trimer. The acceptor LHCII sites behave 
as a sink with a trapping rate that depends on how fast the exciton localized on those "acceptor" sites diffuses to 
a neighboring LHCII trimer or to some other part of the overall photocomplex. In order to study this diffusion 
process, we first focus on a single monomer, initialized in the monomer Hamiltonian ground state \E\) (localized on 
a610, a611, a612). In order to estimate the delocalization time we plot the populations of the various bands in fig. [17] 
(left plot). We use the previously introduced delocalization measure T>(t) which is plotted in the right part of fig. [17] 
The simulations employ a fixed value of the dephasing rate 7^ = \2ps~ 1 that approximately corresponds to ambient 
temperature. The population of the acceptor sites decreases and becomes of the same order of the other populations 
in about lOps. If we take T>(t) = yO + Al expt/ti + A2 expt/t2 to estimate the delocalization process over the whole 
monomer (T>(t) is a site based measure of delocalization) we see that the relevant time scale to, w 2ps (while for the 
initial fast growth we have t\ 350/s). 

The value of t% is of the same order of the inverse of the trapping rate which is characteristic of the FMO k^r ap = lps 
that we used for the monomer in the previous section, and that we will use for the dimer, the trimer and the tetramer 
in the following section. 
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EFFICIENCY: A COMPARISON BETWEEN CLUSTERS OF MONOMERS 



In this section we examine the efficiency and typical transfer time of the single monomer and of groups of 2, 3, 
and 4 monomers. The clusters are build by connecting monomeric subunits via a site-site interaction: as pointed out 
in [l(| there is evidence of a relatively strong coupling in the cluster {6601,6608,6609}, where 6601 belongs to one 



monomer and {6608, 6609} to another one. In [17|, where the LHCII of the spinach is analyzed, the link between the 
6601 and the 6609 sites is estimated to be 42 cmT 1 . We therefore write the overall Hamiltonian of the complex, say 
the trimer, as: 

tt _ \ K rji I rjint . Tjint . rjint / -i q\ 

n trimer — ^ 11 monomer T n 12 T -"23 ' -"31 ' 

1=1,2,3 

where the terms = V660ii.i>609j ( |6601i)(fo609j | + 1 6609^X6601, | account for the Coulomb interaction between the 
Chi 6601 on the i — th monomer and the Chi 6609 on the j — th monomer. Vb6oii,6609j is chosen to be equal to A2cmr 1 
ilOj. We will also consider the case when the coupling 6601 — 6608 is present and equal to 42cm~ 1 . 

The recombination parameter is fixed at V = O.OOlps -1 , which again is the value used for the FMO in Trapping 
is supposed to be similar to the monomeric case; the output sites are now the group of sites a610, a611, 612 on each 
monomer, and the trapping Hamiltonian is J^. fct rap (|a610i)(a610i + |a611j)(a611j| + a612.;){a612i|), where i labels 
the monomers, with ktrap = lps^ 1 . 

The main result of these simulations is that there is a clear evidence for a dephasing- assisted mechanism that 
enhances the transport efficiency of the systems. This mechanism is an example of environmentally assisted quantum 
transport (ENAQT) [5j-|7|. The result is independent of the structure analyzed. Efficiencies and the average transfer 
times have their optimal values in correspondence of a dephasing rate 7^ > 10ps -1 , i.e. the rate corresponding to 
ambient temperature. 



The antenna configuration 



We first describe the efficiency of the various complexes when they act as antennae. We start by focusing on the 
relevant figures of merit for the single monomeric unit. In fig. [TH]we show the efficiency (left) and the average transfer 
time (right) for two different initial states: I-E14}, mostly localized on the stromal side and I-E13), mostly localized 
on the lumenal side. The plots show that the differences highlighted by our analysis of the dynamics in the previous 
section have relevant effects also in terms of the transport efficiency of the monomer. The stromal-stromal pathway 
that starts with I-E14) results in general in a better efficiency and a smaller r over the whole range of dephasing rates. 
In particular, at 7^ > 10 r is ~ 1.5 times smaller than the corresponding lumenal-stromal value. The mismatch in 
the characteristic time of the two pathways arises from the bottleneck present in the lumenal pathway. 
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FIG. 18: Efficiency *?(7^>) (left) and Average transfer time r(7^) (right) for the monomer with different initial states: \E1_3) 
(mostly localized on the lumenal side) and \Ei4) (mostly localized on the stromal side). 

We now describe in detail the results of our simulations for more complex structures. We first examine the case in 
which there is an active sink attached to each of the monomers of a given structure. The initial state used for the 
various simulations is the highest excited state of the Hamiltonian, which is basically a copy over different monomers 
of the eigenstate I-E14) of one single monomer: it is thus mainly localized on the b-stromal Chls of each structure. In 
Fig. [TO] (left) we compare the efficiency of the energy transport for different geometries in presence of only one inter- 
monomeric coupling. While the dephasing assisted mechanism is always present, we see that for physically relevant 
values of the dephasing parameter the change of geometry does not provide significant modifications in the efficiency. 
Fig. HU (right) shows the result of the simulations for t, the characteristic transfer time. The value corresponding to 
the highest efficiency in the transport is around 10-15 ps. This transport time is consistent with that observed for 
the FMO complex, taking into account the fact that each monomer subunit of the LHCII has twice the number of 
chromophores of FMO. 

The results we obtain in this paper are robust with respect to the introduction of static disorder. In Fig |20l we add 
the effects of static disorder (on site-energy values) and compare the efficiency for different numbers of monomers. 
The strength of disorder is taken to be that reported in [26| . The qualitative behavior of the efficiency with disorder 
is the same as that without. The monomer is more affected by disorder than the dimer and trimer, indicating greater 
robustness for the more complicated geometries. This is consistent with the fact that fluctuations due to disorder are 
stronger for systems of smaller size. 

Differences between the various geometries can be observed when the inter-monomeric coupling is supposed to be 
present between the b601 Chi in one monomer and both b608 and b609 in the neighbouring monomer, see fig. [21] 
The efficiency of the clusters of monomers benefits from this kind of coupling: 77 is enhanced with respect the single 
monomer case over a wide range of dephasing values. Moreover, while the dimer is slightly less efficient than the trimer 
and the tetramer structures that have a higher number of traps, our simulations suggest that there is no advantage in 
adding more than three subunits: the trimer behaves just as well as the tetramer. This results could be an indication 
for a functional selection of the trimeric configuration with respect to the other ones: the trimer could be the result 
of an optimization with respect to the "cost" of the structure. 

Interesting differences between the behaviour of the different structures appear when one considers a variable number 
of active sinks. In Fig. [52]we compare the transport in monomeric and dimeric, trimeric and tetrameric antennas with 
only one inter-monomeric coupling, and with different numbers of sinks attached to the available monomeric subunits. 
The initial state in these simulations is always the highest excited state of the structure. The main feature here is 
that the efficiency at physiological temperatures is always greater for the structure with a number of sinks equal to 
the number of monomers. For example, with a single sink attached, the monomer is more efficient than the dimer, 
trimer and tetramer. This behaviour is reasonable since the exciton is initially delocalized over the whole structure 
and in particular over those monomers that do not have any sink attached. A qualitative explanation of the relative 
behavior for the different configurations is based on the competition between an enhancement of the efficiency due to 
a greater number of sinks, and the depletion of efficiency due to a greater number of chromophores, where the exciton 
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FIG. 19: Efficiency r\ (left) and average transfer time r (ps) (right)for the different topologies (monomer, dimer, trimer and 
complex with 4 monomers) with single intra-monomeric coupling for various values of the dephasing rate Initial state: 
highest energy eigenstate of the given structure. 
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FIG. 20: Efficiency r\ for different topologies (monomer, dimer, trimer) in the presence of static disorder, varying the dephasing 
rate 70 . Solid lines represent the averaged value, while dashed lines represent the result for the non-disordered Hamiltonian 
that we use in this work. 

can delocalize and dissipate in the bath, causing a slow down in the funneling process. As already noticed, Fig. H~9l 
the trimer saturates the enhancement of efficiency due to the number of exits: the tetrameric complex with four sinks 
and the trimer complex with three sinks have the same efficiency. 

The wire configuration 

We now analyze the multi-monomer structures when they are used as "wires" . As noted above, the LHCII complexes 
can function as connecting structures between different units of the PSII complex. In the following we analyze this 
case by fixing as the initial state of the dynamics the ground state \E±) of one single monomer in the cluster and by 
connecting a sink to each of the other monomers in the complex. 

Fig. [53] shows the efficiency and the characteristic time of the dimer, the trimer and of a cluster of four monomers 
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FIG. 21: Efficiency rj (left) and average transfer time r (ps) (right)for the different topologies (monomer, dimer, trimer and 
complex with 4 monomers) with two intra- monomeric coupling for various values of the dephasing rate 7^. Initial state: highest 
energy eigenstate of the given structure. 




FIG. 22: Efficiency of the transport for dimer, trimer and tetramer clusters, with one inter-monomeric coupling, acting as 
ANTENNAS In the legend Q23 means that we consider a tetrameric complex where the exciton is captured in the first 
monomer, and can exit only from the II and III monomers. Analoguosly for the trimer (T) and the dimer (D). 



coupled together with only one inter-monomer connection. The simulation suggests again a special role played by 
the trimer. It is significantly more efficient than the dimer, but almost indistinguishable from the cluster of four 
monomers. In the wire configuration, the exciton passing through the dimer has only one trapping site that it can 
reach, while in the trimer, tetramer, and higher order polymers, the exciton entering at the lowest energy site has two 
trapping sites that it can reach - those on the two monomers adjacent to the entrance site. Consequently, the efficiency 
of trapping is higher for the trimer and tetramer than for the dimer: two traps are better than one. Moreover, the 
dynamics for passing from the entrance site to those two adjacent traps are identical for the trimer, tetramer, and 
higher order polymers. Consequently, the efficiency profiles for the trimer, tetramer, and higher order polymers are 
identical. 

We finally describe how the number of attached sinks modifies the overall efficiency of the different quantum wires. 
In Fig. [24]the simulations refer to a situation where the initial state of the structure is the ground-state \Ei) localized 
in the a610, a611, a612 Chls of a single monomer of the structure. Here again we see a behaviour similar to the 
antenna confirguration case: at fixed number of sinks attached the efficiency is higher for those structures with a 
smaller number of chromophores . The trimer is again a limiting case, it is always more or as efficient as the tetramer, 
and, as already shown in fig. (|23p the trimer with two sinks has the same performance of the four-monomer structure 
with three sinks. An analogue comparison holds true for the dimeric vs trimeric cluster. When only one sink is active, 
at the relevant values of dephasing - i.e. where the efficiency is maximal - the dimer can perform slightly better than 
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FIG. 23: Efficiency rj (left) and average transfer time r (ps) (right) for the different topologies (dimer, trimer and complex with 
4 monomers) with single intra-monomeric coupling for various values of the dephasing rate 7^. Initial state \Ei), i.e. ground 
state of a single monomer, localized on the "acceptor" sites a610, a611, a612 . 
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FIG. 24: Efficiency of the transport for dimer, trimer and tetramer clusters, with one inter-monomeric coupling, acting as 
WIRES. In the legend Q23 means that we consider a tetrameric complex where the exciton is captured in the first monomer, 
and can exit only from the II and III monomers. Analoguosly for the trimer (T) and the dimer (D). 

the trimer. From our simulations it is also evident that, within our assumptions for the inter-monomer links, the 
various structures show a directionality of transport. For example the efficiency of the dimer changes significantly 
over a wide range of dephasing values depending on the direction of the energy flow. 

CONCLUSION 

We have analyzed excitonic transport in the primary component of the photosynthetic apparatus in green plants: 
LHCII. By means of a simple decoherence model (Haken-Strobl) , the analysis shows that an exciton initially localized 
on a single chromophore moves through the LCHII photocomplex in a two-step process which clearly signaled by a 
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quantum information motivated measure of derealization. First, over the timescale of a few hundreds of picosec- 
ond, the exciton spreads coherently to neighboring chromophores. This coherent spreading exhibits rapid quantum 
oscillations and entanglement. Second, as the environment decoheres the exciton's position, the exciton diffuses semi- 
cohcrently throughout the complex over a timescale of ten to twenty picoseconds. Although the Haken-Strobl model 
does not include relaxation, we expect this two-step, coherent-semicoherent model to hold for more detailed models 
of the dynamics as well, for example, in the full hierarchy approach of 25| . 

Detailed comparison of different measures of correlations also allows us to identify how the two main downhill 
relaxation pathways unveiled by the experiments (stromal-stromal, lumenal-stromal) can be dynamically described in 
terms of correlations pathways: the correlations among subsy terns of pigments are mostly established on the stromal 
side of the monomeric complex even when the excitation is initially localized on the lumenal side. This behaviour has 
important consequences when the efficiency of the transport is considered: the stromal pathway is in general more 
efficient than the lumenal one. 

In general, the analysis shows that even in the absence of relaxation, pure dephasing induces effective transport 
down the LHCII energy funnel. The transport is efficient and robust in the presence of static disorder, and exhibits 
the characteristic signature of characteristic environmentally assisted quantum transport (ENAQT) - low efficiency at 
low temperature due to transient localization, followed by a robust maximum efficiency at physiological temperature, 
with a falling off of efficiency at very high temperature. In the second part of the paper, we compared the efficiency 
of transport through LHCII structures with different topologies (monomers, dimers, trimers, and tetramers) and 
different configurations (antenna and wire). We find that the efficiency of transport increases as the number of 
subunits increases, saturating at the level of the trimer. These results provide evidence for the functional selection of 
the trimcric configuration. 
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